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ABSTRACT 



The production of the radioactive isotopes '^'^Ti and '^°Co in aU types of supernovae 
is examined and compared to observational constraints including Galactic 7-ray surveys, 
measurements of the diffuse 511 keV radiation, 7-ray observations of Cas A, the late 
time light curve of SN 1987A, and isotopic anomalies found in silicon carbide grains in 
meteorites. The (revised) line flux from '^'^Ti decay in the Cas A supernova remnant 
reported by COMPTEL on the Compton Gamma-Ray Observatory is near the upper 
bound expected from our models. The necessary concurrent ejection of ®'^Ni would also 
imply that Cas A was a brighter supernova than previously thought unless extinction in 
the intervening matter was very large. Thus, if conflrmed, the reported amount of '^'^Ti 
in Cas A provides very interesting constraints on both the supernova environment and its 
mechanism. The abundances of '^'^Ti and '^°Co ejected by Type II supernovae are such that 
gamma-radiation from '^'^Ti decay SN 1987A could be detected by a future generation of 
gamma-ray telescopes and that the decay of '^°Co might provide an interesting contribution 
to the late time light curve of SN 1987A and other Type II supernovae. To produce the solar 
■^■^Ca abundance and satisfy all the observational constraints, nature may prefer at least 
the occasional explosion of sub-Chandrasekhar mass white dwarfs as Type la supernovae. 
Depending on the escape fraction of positrons due to ®'^Co made in all kinds of Type la 
supernovae, a signiflcant fraction of the steady state diffuse 511 keV emission may arise 
from the annihilation of positrons produced during the decay of '^'^Ti to '^'^Ca. The Ca 
and Ti isotopic anomalies in pre-solar grains conflrm the production of '^'^Ti in supernovae 
and that extensive mixing between zones has occurred, but a quantitative model for this 
mixing is presently lacking. 



Subject Headings: gamma rays: theory - nuclear reactions, nucleosynthesis, abundances - 
supernovae: general 



2 



1. Introduction 



One byproduct of Galactic chemical evolution is the production of traces of short to 
medium lived radioactivities having lifetimes ranging from several years to several mil- 
lion years. Detecting the 7-radiation from the decay of these species, either in individual 
supernova remnants or through the glow of accumulated nucleosynthetic products in the in- 
terstellar medium provides a direct way to study the nature of supernovae and to calibrate 
nucleosynthesis theory (e.g., Clayton, Colgate, & Fishman 1969; Clayton 1974, Ramaty 
& Lingenfelter, 1977). In a previous paper we analyzed the production of two long lived 
isotopes ^'^Al and '^°Fe (Timmes et al. 1995b). Here we concentrate on two shorter lived 
isotopes, ^^Ti (half-life of 58 ± 10 yr; Meifiner et al. 1995) and «°Co (half-life of 5.27 
years; Tuli 1990). In addition to their gamma-radiation, these species have additional 
observational consequences for the late time (one year to several centuries) light curves of 
supernovae and for isotopic patterns measured in presolar grains. Moreover, since their 
decay lifetimes are comparable to or shorter than the average time interval between su- 
pernovae, their abundances can serve as important tracers of the most recent supernovae 
in our Galaxy. In this paper we summarize the production levels of these two isotopes in 
all types of supernovae and discuss the observational implications. 

Our starting point is a grid of massive star models recently followed through their 
presupernova evolution by Weaver & Woosley (1996) and exploded, in a parametric way, 
by Woosley & Weaver (1995). Each star was exploded using a piston to give a specified final 
kinetic energy of the ejecta at infinity (typically 1.2 X 10®^ erg). In both the presupernova 
and exploded models the nucleosynthesis of 200 isotopes was determined. These give ejecta 
whose composition, when integrated into a Galactic chemical evolution model, compares 
favorably with what is observed in the Sun and other nearby stars (Timmes, Woosley & 
Weaver 1995a). We address the effects of mass loss by also including the stellar models of 
Woosley, Langer, & Weaver (1993, 1995), which explicitly followed the mass loss that might 
be appropriate for either single or mass exchanging binary stars evolved to supernovae. 
The mass of '^'^Ti and '^°Co ejected is only affected if these stars lose not only all of their 
hydrogen envelopes, but end up with significantly smaller helium cores. Such a reduction 
in helium core mass seems necessary to produce the light curves of Type lb supernovae 
(Ensman & Woosley 1988). 

We also examine the synthesis of '^'^Ti and '^°Co by various models for Type la su- 
pernova. Thermonuclear explosions of accreting carbon-oxygen white dwarfs that have 
approached the Chandrasekhar limit eject a mass of '^'^Ti that is typically, though not al- 
ways, substantially smaller than the yields from Type Il/Ib supernova. Models of helium 
detonations on sub-Chandrasekhar mass carbon-oxygen white dwarfs can synthesize large 
amounts of '^'^Ti (Woosley, Taam & Weaver 1986; Livne & Glasner 1991; Woosley & Weaver 
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1994; Livne & Arnett 1995). If explosion of these types of Type la supernova are not rare, 
they may be the main origin of ■*'*Ca in the Galaxy. 

The main result of this paper is to show consistency between the masses of '^'^Ti and 
'^°Co ejected in current supernova models and several independent sets of observations 
- namely Galactic 7-ray surveys, measurements of the diffuse 511 keV radiation, 7-ray 
observations of Gas A, the late time light curve of SN 1987A, and the isotopic anomalies 
found in silicon carbide grains in meteorites. 

2. The Production of '^'^Ti and '^°Co in Supernovae 

In massive stars, '^'^Ga is produced chiefly, though not exclusively, as the radioactive 
progenitor '^'^Ti made in the "a-rich freeze-out". This process occurs when material ini- 
tially in nuclear statistical equilibrium (NSE) at a relatively low density is cooled rapidly 
enough that the free a-particles present in the NSE ensemble do not have time to merge 
back into the iron group by the relatively inefficient triple-a reaction. Thus, the dis- 
tribution of nuclei cools down in the presence of an anomalously large concentration of 
a-particles (Woosley, Arnett & Glayton 1973; Thielemann, Hashimoto & Nomoto 1990; 
Woosley & Hoffman 1991). The isotope '^°Co is made by neutron capture in the neon and 
oxygen layers, so its synthesis is less sensitive to uncertainties in the explosion model. 

Yields of '^'^Ti and '^°Co from the (solar metallicity) Type II supernova models of 
Woosley & Weaver (1995) are shown in Fig. la, and from the Type lb models of Woosley, 
Langer & Weaver (1995) in Fig. lb. Production in Type lb supernovae is more uni- 
form because the models converge to a common presupernova mass in the narrow range 
2.26 - 3.55 Mq. All of the models shown in Figs, la and lb, whose ejecta all have nearly 
]^q5i gj,g q£ J^inetic energy at inflnity, predict '^'^Ti yield between 10~® and 1.5 X 10"'* Mq, 
with typical values ~ 3 X 10~® Mq for the Type II models and ~ 6 X 10~® Mq for the 
Type lb models. Note that the flnal mass of the helium core is a more appropriate pa- 
rameterization than either the main-sequence or initial helium core mass of the star. The 
supernova models also predict a yield of '^°Co between 10"'^ to 10"'* Mq, with typical values 
being ~ 3 X 10~® Mq for the Type II models and ~ 2 X 10"'^ Mq for the Type lb models. 

Producing more than 10"'* Mq of '^'^Ti is difficult in core collapse models, but making 
much less is easy. Indeed, ejection of any '^'^Ti is sensitive to how much mass falls back 
onto the remnant (Woosley & Weaver 1995). The mass of '^'^Ti external to the piston and 
the mass ejected is shown in Table 1 as a function of the main-sequence progenitor mass. 
The difference is the portion that falls back onto the neutron star. Entries in the table are 
for the models where the ejecta's kinetic energy is a nearly constant 10®^ erg. The masses 
external to the piston are probably realistic upper bounds. In all cases, the mass external 
to the piston is greater than or equal to the ejected mass, with equality holding only for 
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the least massive stars (11 and 12 Mq in Table 1). To achieve a nearly constant kinetic 
energy of the ejecta, the explosion energy in the M > 25 Mq presupernova models must 
be steadily increased in order to overcome the increased binding energy of the mantle. In 
the most massive stars (35 and 40 Mq in Table 1) the binding energy and density profile 
causes nearly all the '^'^Ti produced to fall back onto the compact remnant (see Woosley 
& Weaver 1995 for a discussion of this fall back). Unless the explosion mechanism, for 
unknown reasons, provides a much larger characteristic energy in more massive stars, it 
appears likely that stars larger than about 30 Mq will have dramatically reduced '^'^Ti 
yields and leave massive remnants (M > 10 Mq) which become black holes. If, however, 
the explosion is energetic enough (perhaps overly energetic) more '^'^Ti is ejected. Clearly, 
other than the upper bounds of Table 1, the '^'^Ti synthesis in models of very massive stars 
is quite uncertain. 

In the solar metallicity 25 Mq model, the isotope '^°Co, like '^°Fe, is made chiefly in the 
oxygen and neon shells in the last stages of presupernova evolution (see Fig. 2a) and to a 
lesser extent in a thin layer at the base of the He-burning shell. The production of '^°Fe by 
neutron capture in the explosive helium burning shell is larger than the production of '^°Co 
(made also by neutron capture, see Fig. 2b), since the ^'^'^®Fe seed nuclei are more abundant 
than the ®®Co seed nuclei. The difference between production by neutron exposure in the 
helium shell and in the convective neon-oxygen shell may be shown by examining the 
integrated totals. Exterior to 6.53 Mq (the bottom of the He shell) there is 1.9 X 10~^ Mq 
of «°Fe and 3.4 x 10"^ Mq of «°Co in the 25 Mq presupernova model and 1.2 X 10"^ Mq 
of ®°Fe and 4.5 X 10"® Mq of ®°Co after the explosion. Between 2.07 Mq (the future mass 
cut) and 6.53 Mq the presupernova has 7.8 X 10"® Mq of ®°Fe and 4.7 x 10"^ Mq of 
®°Co, while after the explosion there is 9.4 X 10"® Mq of ®°Fe and 4.9 X 10"^ Mq of ®°Co. 
Thus, the total ejected mass of ®°Fe is 2.1 X 10"^ Mq and 5.4 X 10"^ Mq of ®°Co as 
given in the tables of Woosley & Weaver (1995). In addition to its sensitivity to explosion 
characteristics and the strength of the s-process, the amount of '^°Co made here might be 
altered by a more careful treatment of convection. Bazan & Arnett (1994) flnd that plume 
structures dominate the two dimensional velocity fleld of convective oxygen shell burning 
and that signiflcant local hot spots of nuclear burning exist due to mixing beyond the 
boundaries deflned by mixing-length theory. The composition inhomogeneities and local 
burning are likely to change the quantitative yields of several isotopes from a single star 
as compared to one-dimensional calculations. Yields from a well deflned progenitor mass 
may even be mildly chaotic. 

Since we are dealing with neutron rich nuclei near the boundary of the network em- 
ployed by Woosley & Weaver (1995), it was desirable to check the results using a larger set 
of reactions. The solar metallicity 15 and 25 Mq stars were evolved with a larger nuclear 
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reaction network (500 isotopes, up to and including the isotopes of ruthenium). Also in- 
cluded was a temperature dependent beta decay rate for '^°Co which effectively decreases 
its synthesis during helium burning. After the entire evolution (presupernova + explosion) 
the mass of '^'^Ti ejected by the 15 Mq star decreased 4.9% and the mass of '^°Co ejected 
decreased by 25%. For the 25 Mq model, '^'^Ti increased by 25% while '^°Co decreased by 
-3.1%. We conclude that the error in the abundances of these isotopes due to the smaller 
size of network used by Woosley & Weaver (1995) is negligible. 

Deflagration models for Type la supernovae in which carbon ignites at or near the 
center of a carbon-oxygen white dwarf of 1.39 Mq produce 1.8 X 10~® Mq of '^'^Ti and 
2 X 10"® Mq of «°Co (Nomoto, Thielemann & Yokoi 1984; Thielemann, Hashimoto & 
Nomoto 1986). While the synthesis of '^°Co is negligible, the yield of '^'^Ti in the W7 Type 
la model is within a factor of two of the Type Il/Ib models (see Fig. 1). However, core 
collapse supernovae occur about 6 times more frequently than thermonuclear supernovae 
in the Galaxy (van den Bergh & Tammann 1991; van den Bergh & McClure 1994). If 
all Type la events eject 1.8 X 10~® Mq of '^'^Ti and the average Type H/Ib event ejects 
3 X 10~® Mq, then standard deflagration Type la supernovae account for only ~ 8% of 
the Galactic '^'^Ca abundance. Some other deflagration models (Woosley & Eastman 1995) 
produce even less '^'^Ti, while still accounting for the majority of the spectral and light 
curve observations. It appears that standard Type la deflagration models do not make an 
appreciable fraction of the Galactic '^'^Ca abundance. 

Several recent stellar evolution studies suggest an alternative progenitor for Type la 
supernovae, a 0.6 - 0.9 Mq carbon-oxygen white dwarf that merges by gravitational ra- 
diation with a helium main-sequence star and accretes helium at a rate of several times 
10~® Mq yr~^ (Iben & Tutukov 1991; Limongi & Tornambe, 1991; Munari & Renzini, 
1992; Tutukov, Yungelson & Iben 1992; Iben & Livio 1993; Tutukov & Yungelson 1995). 
When 0.15 - 0.20 Mq of helium has been accreted in this model, a helium detonation is 
initiated at the base of the accreted layer. This helium detonation compresses the carbon- 
oxygen core and triggers a detonation there as well (Nomoto 1980, 1982; Woosley, Weaver 
& Taam 1980; Woosley, Taam & Weaver 1986; Livne 1990; Livne & Glasner 1991; Woosley 
& Weaver 1994; Woosley & Eastman 1995). The iron group nucleosynthesis is acceptable in 
these models (0.3 - 0.9 Mq of ^^Fe ejected), with 2.9 x 10"^ - 3.9 x 10"^ Mq of ^^Ti ejected. 
Production factors for '^'^Ca relative to its solar mass fraction thus range from 200 - 3000. 
The isotopes 46,47,48rp- siy ^^sq have large production factors in these models. The 
synthesis of '^'^Ti in these models is by a high temperature variety of incomplete explosive 
helium burning. Since the material never achieves NSE, it is not technically an alpha-rich 
freeze out. Depending upon how frequently these types of sub-Chandrasekhar mass white 
dwarfs explode, the large production factors suggest that these types of thermonuclear 
events might be the principal origin of '^'^Ca in nature. Sub-Chandrasekhar models also 
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make significant amounts of '^^Ti and ®^V, which are also deficient in the steUar-chemical 
evolution models. 

3. Detection of '^'^Ti 

The radioactive isotope '^'^Ti decays to '^'^Sc with an accompanying emission of 7-rays at 
67.85 and 78.38 keV. The most recent measurement of the '^'^Ti half-life gives a preliminary 
value of 58 years with a 1 sigma error of 10 years (Meifiner et al. 1995). These values 
cover the range (46.4 to 66.6 yr) of previously reported values (Wing et al. 1965; Moreland 
& Heymann 1965; Frekers et al. 1983; Alburger & Harbottle 1990). The sources of the 
large uncertainty in the ongoing Meifiner et al. (1995) measurement appear to be well 
understood and the uncertainty should drop in the near future to well below 10%. The 
isotope ■^■^Sc then decays to the stable isotope ■*'*Ca with a half-life of 5.7 hr, emitting a 
1.157 MeV 7-ray. This line has been a recognized target for 7-ray astronomy for quite 
some time (e.g., Clayton et al. 1969) 

Direct detection of '^'^Ti is possible through 7-ray (or hard x-ray) observations of in- 
dividual supernova remnants or small groups of extinguishing OB associations. Additional 
means for inferring the abundance of '^'^Ti produced in supernovae are the analysis of late 
time supernova light curves and graphite grains from meteorites. In the remainder of 
this section we examine the evidence from each of these lines of inquiry in view of the 
production levels described in §2. 

3.1 ■^■^Ti in the Galaxy 

Since the half-life of '^'^Ti is comparable to the average time between Galactic super- 
novae, its 7-ray emission should appear as a small collection of Galactic point sources. 
Of the six known local supernovae (which probably comprises a complete sample within 
a distance of ~ 3 kpc) only SN 1006 is considered to have been a Type la event. Gas 
A and Tycho were most likely Type lb supernova, while the Grab, SN 1181 and Kepler 
were probably Type II supernovae. Of these six, only the '^'^Ti abundance in Gas A has 
been pursued aggressively and is discussed in §3.2. Evidence for Galactic '^'^Ti emission 
was sought by the large field of view spectrometers aboard the High Energy Astronomy 
Observatory (HEAO) and Solar Maximum Mission (SMM) satellites. Neither search for 
young, unknown and visually obscured supernovae remnants produced a positive '^'^Ti de- 
tection. Mahoney et al. (1992) estimated a 99% confidence level upper limit of 2 X lO"'* 
photons cm~^ s~^ for the 67.85 and 78.38 keV line fluxes using the HEAO survey. Leising 
& Share (1994) placed a 99% confldence level upper limit of 8 X 10~® photons cm~^ s~^ 
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for the 1.157 MeV line flux in the general direction of the Galactic center by analyzing 10 
years of SMM data. 

Detailed stellar-chemical evolution models of the Galaxy (e.g., Timmes et al. 1995a) 
produce only about 1/3 of the solar '^'^Ca abundance, with most of the production at- 
tributable to Type II supernova. This model provides satisfactorily reproduces the solar 
abundances for the vast majority of isotopes lighter than gallium, dwarf star elemental 
abundance histories, estimates of the present average Galactic supernova rates, spatial 
distribution and flux levels expected from ^'^Al and '^°Fe (Hoffman et al. 1995; Timmes 
et al. 1995b) and the present interstellar medium ^'^Al/^^Al ratio. Note that the sum of 
■^■^Ca made as itself and '^'^Ti synthesized were used in the calculation of the solar '^'^Ca 
abundance. The mass of '^'^Ca produced as itself is comparable to the '^'^Ca synthesized as 
radioactive ■*'*Ti. 

To explain the solar '^'^Ca abundance, either the Type II '^'^Ti yields from Woosley & 
Weaver (1995) are too small (e.g., because of too much fall back in the more massive 
stars), the calculated supernova rates are too small, or there is another source of '^'^Ca in 
the Galaxy. Since '^'^Ti is produced in the innermost ejecta of the supernova, it is strongly 
dependent on the explosion details. This delicate dependence can be seen in Fig. 1 with 
the somewhat irregular behavior on progenitor mass. Explosions that one may regard as 
overly energetic, and neglecting any fallback material, alleviates most of the factor of 3 
discrepancy (see Table 1). The iron yields in this case are, of course, much larger and may 
pose a problem for duplicating the abundance histories observed in dwarf stars. 

Thielemann, Nomoto & Hashimoto (1996) have also examined the detailed nucleosyn- 
thesis of core collapse supernovae, and they flnd larger amounts of '^'^Ti is ejected. Their 
13 Mq star ejects 2.4 - 5.0 x 10"^ Mq of ^^Ti, a 15 Mq star 7.2 - 12.0 x 10"^ Mq, a 20 Mq 
star 17.0 x 10~® Mq, and a 25 Mq star 2.1 x 10~® Mq. These values are sometimes similar, 
sometimes not, to those given in Table 1. Some of the differences between the Woosley & 
Weaver and Thielemann et al. yields are caused differences in the nuclear reaction rates. 
Hoffman et al. (1996) show that use of the Thielemann et al. reaction rates in the Woosley 
& Weaver stellar models produce nearly the same yields for the vast majority of isotopes. 
For the speciflc case of '^'^Ca, Hoffman et al. (1996) show a factor of 2 spread between the 
two groups due to differences in the relevant reaction rates. 

Another reason for the differences in the yields is tied to the different adiabatic paths 
followed in the explosion, since the progenitor structure is quite similar between the two 
groups (Hoffman et al 1996). Thielemann et al. explode their stellar models by depositing 
thermal energy deep in the neutronized core (which later becomes the neutron star). The 
energy deposited gives a larger entropy to the innermost zones than what a piston driven 
explosion would impart and, in principle, a thermal explosion would eject more material. 
A larger entropy ensures a more vigorous alpha-rich freeze out and, thus, a larger '^'^Ti 
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production. However, Thielemann et al. sum the ejecta from the outside in and place the 
mass cut (which is artificial in all models) at the position where sufficient ®'^Ni is produced 
to explain the observations. Woosley & Weaver, on the other hand, place the piston at 
a suitable mass cut and then follow the explosion. It is the difference between injecting 
momentum or energy in modeling the explosion that leads to the two groups following 
different adiabatic paths (also see Aufderheide, Baron & Thielemann 1991). What real 
self-consistent explosion models do, exploding via neutrino heating and multidimensional 
convection, might be reffective of one dimensional momentum driven explosions, one di- 
mensional energy driven explosions, or some intermediate case. Hence, the differences 
between the two groups in the amount of '^'^Ti ejected is an example of the spread one 
obtains due to uncertainties in modeling the Type H explosion mechanism. 

A leading candidate for an additional source of '^'^Ti would be the sub-Chandrasekhar 
mass white dwarf models for Type la supernovae discussed in the previous section or other 
Type la models where a detonation wave passes through a surface layer of helium. Nature 
may also prefer sub-Chandrasekhar mass Type la supernovae to make the isotopes '^^Ti and 
®^V, which are also deficient by a factor of 3 in the stellar-chemical evolution models. A 
larger '^'^Ti abundance from any of these three possibilities (larger yields, larger supernova 
rate, additional source) could account for a significant fraction of the strong diffuse Galactic 
511 keV emission due to electron-positron annihilation (Bussard, Ramaty, & Drachman 
1979; Woosley & Pinto 1988; Lingenfeffer, Chan, & Ramaty 1993). 



3.1.1 The diffuse Galactic annihilation radiation 



Observations of the 511 keV line from positron annihilation show a steady diffuse 
component from the Galactic disk superimposed upon a time-variable point source located 
near the Galactic center (Johnson, Harden & Haymes 1972; Leventhal, MacCallum & 
Stang 1978; Niel et al. 1990; Gehrels et al. 1991; Skibo, Ramaty, & Leventhal 1992; 
Leventhal et al. 1993; Purcell et al. 1993; Ramaty, Skibo. & Lingenfelter 1994). Decay 
products from radioactive decay of ^'^Al, '^'^Ti and ®'^Co produced by supernovae, novae or 
the winds of Wolf-Rayet stars have been proposed as a source for the diffuse component 
(Ramaty & Lingenfelter 1979; Woosley & Pinto 1988; Lingenfelter et al. 1993), along 
with smaller contributions from the decay of pions produced by cosmic ray interactions 
(Ramaty, Kozlovsky & Lingenfelter 1979; Ramaty & Lingenfelter 1991), energetic solar 
fiares (Mandzhavidze & Ramaty 1992), pulsars (Sturrock 1971) and gamma-ray bursts 
(Lingenfelter & Hueter 1984). However, none of the previous estimates for the contribution 
from radioactive nuclei produced in supernova enjoy the number and/or quality of the 
stellar models employed here. 
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Two positrons are produced in 95% of the decays from '^'^Ti to ■*'*Ca, while 85% of the 
decays from ^'^Al to ^'^Mg produce a single positron. That portion of pair annihilation 
radiation due to the decay of '^'^Ti and ^'^Al balance the steady state production of these 
isotopes. Once the positrons are in the interstellar medium, they have a mean life against 
annihilation of at least 10® yr (Bussard et al. 1979; Lingenfelter et al. 1993; Chan & 
Lingenfelter 1994; Ramaty et al. 1994). This is much larger than the average time interval 
between supernovae, so the steady state assumption is valid. 

For the exponential disk + 1/r^ bulge Galactic mass distribution used by Timmes et 
al. (1995a) and requiring a full production of the solar abundance of ■*'*Ca, we calculate 
annihilation rates of 2.9 X lO'*^ positrons s~^ due to ^'^Al, 4.5 X lO'*^ positrons s~^ due to 
■^■^Ti, and 1.6 X lO'*^ positrons s~^ from ®'^Co. The total annihilation rate is then 2.3 X lO'*^ 
positrons s~^ with 13 % of the 511 keV emission explained by the decay of ^'^Al in the 
interstellar medium. Lingenfelter et al. (1996), using a similar mass distribution, estimated 
a total annihilation rate of 1.5 it 0.5 X lO'*^ positrons s~\ while measurements of the diffuse 
1809 keV ^'^Al 7-ray line emission constrain the positron contribution due to ^'^Al to 16 it 
5% (Mahoney et al. 1982, 1984; Share et al. 1985; Gehrels et al. 1993; Skibo et al 1992; 
Diehl et al. 1995) 

In this simple picture, the '^'^Ti yields from Type la and Type Il/Ib were artificially in- 
creased by a factor of 3 so that a solar ■*'*Ca abundance was obtained at a time 4.6 Gyr ago 
at a distance 8.5 kpc from the Galactic center (see the discussion in §3.1). The supernova 
rates calculated by the stellar-chemical evolutions were 3.1 per century for Type II -|- lb 
supernova and 0.53 per century for Type la supernova, in good agreement with the esti- 
mates of van den Bergh & Tammann (1991) and van den Bergh & McClure (1994). The 
annihilation rate due to ®'^Co assumed an average escape fraction of 3%. However, ®'^Co 
has a short mean lifetime of 111.4 days and positrons are only produced in 19% of the 
®'^Co decays so that a large fraction could annihilate inside the supernova ejecta before the 
mean free path becomes long enough for the positrons to escape. In addition, a magnetic 
field of 10"'^ gauss gives a gyroradius smaller than the homologous expansion radius at 
111.4 days. Thus, it is possible that the escape fraction of positrons from ®'^Co is an order 
of magnitude smaller than the 3% used above. 

In the limiting case that the escape fraction of positrons from the decay of ^^Co is 
zero, and the '^'^Ti mass ejected by supernovae is pushed to the point of over-producing 
the solar '^'^Ca abundance by a factor of 2, the annihilation rate due to '^'^Ti increases to 
lO'*^ positrons s~^ The total annihilation rate in this case is 1.3 X lO'*^ positrons s~^ with 
23% emission explained by the decay of ^'^Al. The error bars on the total annihilation 
rate and the ^'^Al contribution (unresolved angular distribution, subtraction of the time- 
variable contribution, amount of discrete foreground source contamination) as well as the 
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uncertainties in the supernova models thus aUow the possibihty that a significant fraction, 
if not most, of the smooth 511 keV emission is from the annihilation of '^'^Ti positrons. 

3.2 ^^Ti in Cas A 

The Cas A supernova remnant is relatively close (~ 3 kpc), young (~ 300 yr) and 
wide (physical diameter ~ 4 pc), making it one of the premier sites for studying the 
composition and early behavior of a supernova remnant. It may only be equaled as SN 
1987A unfolds. During July 1992 and February 1993 Cas A was observed (/=111.73°, 
6=-2.13°) with the COMPTEL telescope aboard the Compton Gamma-Ray Observatory 
(CGRO), and showed evidence for line emission at 1.157 MeV from the decay of '^'^Ti at 
a significance level of ~ 4cr (lyudin et al. 1994). A total of 330 photons were attributed 
to the 1.157 MeV feature. Adopting a distance of 2.8 kpc, the initially measured line flux 
of 7.0 lb 1.7 X 10~® photons cm~^ s~^ (detection limit of 10~® photons cm~^ s~^) converts 
into a ■^■^Ti mass ejected during the Cas A supernova explosion of 1.4 - 3.2 X lO"'* Mq, 
depending on the precise date of the event and the mean lifetime of '^'^Ti (lyudin et al. 
1994). This was the flrst (and so far only) time a positive detection of 7-ray emission from 
the decay of '^'^Ti has been reported, either in the debris left by individual supernova or 
through point sources of Galactic emission. 

A recent analysis of additional COMPTEL observations of Cas A were unable to at- 
tribute a comparable number of photons to the 1.157 MeV feature, giving a lower flux 
of 4.2 lb 0.8 X 10~® photons cm~^ s~^ (Schonfelder et al. 1995), corresponding to slightly 
smaller ^^Ti mass of 1 - 2 X 10"^ Mq. 

During July 1992 and February 1993 the Galactic plane in the region near Cas A was 
also observed with OSSE CGRO, and showed no flrm evidence of emission from the three 
strongest decay lines (67.85 keV, 78.38 keV and 1.157 MeV) of ^^Ti (The et al. 1995). That 
is, OSSE did not conflrm the COMPTEL detection of the '^'^Ti decay lines. From simulta- 
neous flts to the three 7-ray line fluxes and reasonable estimates of the systematic errors. 
The et al. (1995) derive a 99% confldence upper limit of 5.5 X 10~® photons cm~^ s~^ to 
the flux in each line. Adopting a distance of 2.92 kpc to the Cas A remnant, the upper 
limit line flux corresponds to 1.7 - 5.0 X lO"'* Mq of '^'^Ti ejected during explosion, de- 
pending again on the age of Cas A and the '^'^Ti mean lifetime. Considering the statistical 
uncertainties only, the COMPTEL (lyudin et al. 1994) and OSSE (The et al. 1995) ^^Ti 
flux measurements are formally consistent with each other only at the 1.8% confldence 
level. 

More recently however, OSSE observations of Cas A (The et al. 1996) have tripled 
their exposure time. The best simultaneous flt to the flux of the three '^'^Ti lines yields 
1.7 ± 1.4 X 10-5 photons cm'^ s'^ (The et al. 1996). Thus, the total COMPTEL and 
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OSSE observations are now less in conflict, and may be converging on a lower value than 
COMPTEL originally reported. At the la level the two experiments are in agreement near 
3.5 X 10-5 cm-2 s-\ implying 1.0 x 10"^ Mgof ^^Ti. 

The larger values of '^'^Ti ejected during explosion that are inferred from the COMPTEL 
observations would be difficult to accommodate in the supernova models; the smaller ones 
might not (see Figure 1). Hoffman et al. (1995) have pointed out that a large mass of 
■^■^Ti cannot be ejected unless a large quantity of ®'^Ni is also ejected, and this implies a 
bright supernova. How bright? An abundance of '^'^Ti as large as reported (~ lO"'* Mq) 
would probably mean the ejection of at least 0.05 Mq of ®'^Ni and possibly considerably 
more (Woosley et al. 1995). With or without a hydrogen envelope, the supernova would 
then have a peak luminosity brighter than the lO'*^ erg s~^ observed for SN 1987A. For a 
distance of 2.8 kpc to Cas A, this would imply a peak absolute magnitude of -4, assuming no 
interstellar absorption. Cas A was certainly not detected as such a bright object. Ashworth 
(1980) has argued for a possible eyewitness observation of the Cassiopeia supernova by John 
Flamsteed (flrst astronomer royal of England) in 1680 when the new star had an apparent 
brightness of 6th magnitude. The absence of widespread reports of the event probably 
means that the peak visual magnitude was fainter than 3rd magnitude (van den Bergh 
1991; Da Silva 1993). A peak absolute magnitude of -4 certainly would have been recorded 
by both Eastern and Western observers. 

However, there may be large reddening corrections. Minkowski (1957) and Searle 
(1971) estimate a visual extinction of 4.3 magnitudes in the direction of Cas A, while 
Peimbert & van den Bergh (1971) suggest a visual reddening correction of 5 - 7 magni- 
tudes due to the circumstellar material surrounding the supernova remnant. If Cas A 
was embedded in a dusty region or experienced signiflcant mass loss which condensed into 
dust grains before the explosion, the extinction could have been even larger. Recent mea- 
surements of the X-ray scattering halo around Cas A (Predehl & Schmitt 1995) from the 
Rontgen satellite (ROSAT) offer some compelling evidence for a larger reddening correc- 
tion. The ROSAT measurements show that the interstellar absorption column depth of 
Nh is a factor of ~ 2 larger than the column depth of Nh that normally corresponds to an 
optical extinction of Ay = 5. The material causing the Nh excess appears to be metal- 
rich, dust poor and distributed close (relative to its distance) to the remnant (Predehl 
& Schmitt 1995). If this material had a typical gas/dust ratio at the time of the Cas A 
explosion, then an additional visual extinction of 5 magnitudes may apply. A total visual 
extinction of 10 magnitudes is near the value needed to reconcile the Flamsteed detection 
and our estimate of the peak magnitude based on nucleosynthesis arguments. It is thus 
conceivable that a dusty shell prevented the widespread detection of this supernova, and 
that the supernova shock wave destroyed the dust as it propagated through the debris and 
the material surrounding the Cas A supernova (Hartmann et al. 1996). 
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A progenitor more massive than 10 Mq has been suggested for Cas A by analysis 
of the fast moving optical knots in the supernova remnant (Brecher & Wasserman 1980; 
Ellison et al. 1994), while estimates of the mass emitting X-rays range from 2.4 - 15.0 Mq 
(Fabian et al. 1980; Tsunemi et al. 1986; Jansen et al. 1988; Ellison et al. 1994). The 
three-dimensional structure of the remnant has been explored (Reed et al. 1995). These 
estimates suggests that if Cas A was a Type II event, the progenitor mass was between 
12 and 25 Mq. Combined with the COMPTEL/OSSE observations of ~ 10"^ Mq of ^^Ti 
being ejected, implies that there was little or no fall back during the explosion. The actual 
mass at the time of the explosion is critical. A 20 - 25 Mq supernova might eject ~10~'* Mq 
of ■^■^Ti and still leave a remnant of baryonic mass 2.0 Mq that might become a black hole 
(Woosley & Weaver 1995). A presupernova lighter than 19 Mq would almost certainly 
leave a neutron star (Timmes, Woosley, & Weaver 1996). A Type lb progenitor for Cas 
A (e.g., Fesen, Becker & Blair 1987; van den Bergh 1991; Fesen & Becker 1991) even if 
it had a large main sequence mass, might have lost sufficient mass as a Wolf-Rayet star 
to leave a neutron star remnant. Thus the nature of the compact remnant in Cas A, 
from a theoretical view, is ambiguous. Observationally, X-ray imaging does not reveal the 
synchrotron nebula one might expect around a neutron star (Murray et al. 1979; Pravado 
& Smith 1979; Becker et al. 1980; Jansen et al. 1988; Ellison et al. 1994; Hartmann et al. 
1996). Deep images taken under excellent seeing conditions have set strict magnitude limits 
(M/ > 8.9, M.R ^ 9.3) on the presence of a stellar remnant (Van Den Bergh & Pritchet 
1986). The region around the center of Cas A simply appears void of any detectable stars. 
So observationally, a black hole may have formed. This would make the ejection of a lot 
of ■^■^Ti more difficult, but not impossible. The black hole mass would need to be not too 
far above the critical mass, probably less than 2.0 Mq (gravitational mass). 



3.3 ^^Ti in SN 1987A 



One of the advantages of a nearby supernova like SN 1987A is the opportunity to follow 
the light curve and spectrum to much later times than has been possible in any other 
event. An important effect of '^'^Ti in SN 1987A is to provide a floor to the bolometric light 
curve that is proportional to the uncertain abundance produced in the explosion (Pinto 
& Woosley 1988; Woosley, Pinto, & Hartmann 1989; Kumagai et al. 1989; Arnett et al. 
1989). The thermal luminosity deriving from '^'^Ti decay to '^'^Ca is 

^44 = 4.1 X 10^' [1 - exp (-/€44<^o(to A)') + 1-3] X 

ewi-t/r,,) [y;^^-^,) ergs- , (1) 
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where (po = 7.0 X 10'* g cm ^ is the column depth at some fiducial time to = 10'^ s, and 
K44 = 0.04 cm^ is an estimate of the eff'ective opacity (Pinto & Woosley 1988; Woosley, 
Pinto & Hartmann 1989). 

Observed and theoretical bolometric light curves of SN 1987A are shown in Fig. 3a 
for the 500 to 1500 day time period and in Fig. 3b for the 500 to 3500 day period. In 
each case the observed bolometric luminosity (e.g., Suntzeff et al. 1992) is shown as the 
filled circles. There are no data points past day 2000 since most the emission is in the 
far infrared wavelength region where it is not easily observed. A rather large, but not 
unphysical, ejected mass of M('*'*Ti) = 1.0 X lO"'* Mq was used. Changes in the amount of 
■^■^Ti ejected produce a linear shift of the light curve (see equation 1). The labeled dashed 
lines indicate the 7-ray luminosity of the various radioactive isotopes and the solid line is 
the total luminosity, assuming that radioactive decay is the sole power source. Radioactive 
■^■^Ti tends to dominate prompt radioactive contributions to the bolometric light curve after 
about 1500 days due to the half-life of '^'^Ti and local deposition of the positron kinetic 
energy ('^'^Sc to ■*'*Ca). Figure 4 shows the small effect of the uncertainty in the effective 
opacity K44. Factors of two variation from the value used above are shown as dashed 
lines. The figures suggest that we are now entering an epoch in the supernova's life when 
the dominant energy source, exclusive of input from a pulsar, accreting compact object 
or circumstellar interaction, should be the decay of radioactive '^'^Ti with a small current 
contribution from '^°Co (discussed in §4). This energy input may be difficult to uniquely 
resolve, however, if the atomic processes that convert nuclear decay energy into optical- 
infrared luminosity are no longer operating in steady state (Clayton et al. 1992; Fransson 
& Kozma 1993). 

For a distance of 50 kpc to SN 1987A and a half-life of 58 yr, 5 x 10"^ Mq of ^^Ti 
would produce a 7-ray line fiux of ~ 1 - 2 X 10"'^ photons cm~^ s~^ This line fiux is 
too small for CGRO and probably too small for International Gamma-Ray Astrophysics 
Laboratory (INTEGRAL) instruments, but large enough that it might be detected in the 
next century by post-INTEGRAL experiments. 

3.4 ■^■^Ti in Meteorites 

Interstellar silicon carbide isolated from the Murchison carbonaceous meteorite is 
anomalous in its C, Si, Mg, Ti, Sr, Ba, Nd and noble gas isotopic compositions (Zin- 
ner. Tang k Anders 1989; Ott k Begemann 1990; Lewis, Amari k Anders 1990, 1994; 
Hoppe et al. 1994; see also Anders k Zinner 1993). Most silicon carbide grains carry the 
signature of hydrogen and helium burning with s-process trace elements. A few of the 
grains, however, exhibit extremely exotic isotopic compositions, distinct from the major- 
ity of the silicon carbide grains while being indistinguishable in their morphology. These 



14 



grains were named "type X" grains (Zinner et al. 1991; Amari et al. 1992) and show large 
excesses of ^^C relative to ^^C (up to 76 times solar), relative to ^'^N (up to 22 times 
solar), 26A1/27A1 ratios up to 0.6, depletion in ^'''^°Si relative to ^^Si (up to 77%), and ^^Ti 
excesses up to 95% relative to '^^Ti (Amari et al. 1992; Nittler et al. 1996; Hoppe et al. 
1996; Nittler 1995 private communication). Ion micro-probe mass spectrometer analyses of 
one X type grain from the Murchison carbonaceous meteorite showed a large ■*'*Ca excess 
(relative to '^'^Ca) of 300% (Amari et al. 1992). New (unpublished) data on another 9 
silicon carbide X type grains and a few graphite grains have inferred 44^j^48rp- j^g^^jgg -^yl^ich 
range from 0.0033 to 0.058 (Amari et al. 1995; Nittler et al. 1995), while in another silicon 
carbide grain this ratio is an astonishingly large 0.37 (Hoppe et al. 1996). 

Assuming that grains form in supernovae, the Ca and Ti isotopic anomalies confirm 
that production of ■*'*Ca occurs during explosive nucleosynthesis in supernovae. (Or as the 
meteoritic community prefers, assuming supernovae produce ■*'*Ca the isotopic anomalies 
strongly suggests grain condensation from supernova ejecta.) The requirement that '^'^Ti, 
produced in the deepest layers of the supernova to be ejected, be present in an environ- 
ment that is carbon rich also requires considerable mixing of different zones prior to or 
during the grain formation epoch since carbon exceeds oxygen only in the helium shell. 
Large ■*'*Ca anomalies had been predicted by Clayton (1975) based on the condensation of 
radioactive ■*'*Ti. Clayton envisioned production of '^'^Ti in explosive helium, oxygen, and 
silicon burning. It is synthesized therein, sufficiently for production of '^'^Ca excesses in 
objects having an enhanced Ti/Ca ratio, but those zones do not provide the bulk of '^'^Ti 
synthesized in massive stars. It came as a major surprise to find these titanium anomalies 
in carbide grains, since none of the explosive zones in question contain significant amounts 
of carbon. While it has recently been realized that the ejecta from supernova are subject 
to a variety of macroscopic dynamical instabilities (chiefly Rayleigh- Taylor and Richtmyer- 
Meshkov) that induce mixing between zones (e.g., Herant & Woosley 1994), the need for 
mixing at a microscopic level is still surprising. Until calculations are done that include a 
realistic description of the turbulent cascade to small scales and until our understanding of 
grain chemistry in young supernovae advances further, we are unable to draw quantitative 
conclusions regarding the amount of '^'^Ti required to explain these observations. 

4. Detection of ^°Co 

Each beta decay of '^°Co to '^°Ni (half-life of 5.27 years; Tuli 1990) produces 7-rays of 
1.332 MeV and 1.173 MeV for a total 7-ray energy of 2.505 MeV. The maximum kinetic 
energy of the emitted electron is Q = 0.3178 MeV. We assume that the electron deposits 
its energy locally even when the gamma-rays escape. Lacking a detailed Kurie plot for 
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this nucleus, an estimate of the average kinetic energy of the neutrino is given by the 
approximate formula 

^ ' V 2 / V wy V 4W 9wy ^ ^ 

where 

W=l + -^ = -^. (3) 

The average electron kinetic energy is then the maximum kinetic energy Q minus the 
average neutrino kinetic energy. One concludes that the beta decay electron contributes 
about 7% of the 7-ray energy from '^°Co. Since the electron's kinetic energy is deposited 
locally, however, it remains an energy source in the supernova remnant for several half-lives 
of «°Co. 

At present, detection of '^°Co is limited to Galactic point source 7-ray surveys and the 
analysis of late time supernova light curves. Leising & Share (1994) searched almost 10 yrs 
of data from the SMM Gamma- Ray Spectrometer for evidence of 7-ray line emission from 
the decay of '^°Co. Although Leising & Share found no direct evidence of emission, they 
placed a 99% confidence level upper limit of 8 X 10~® photons cm~^ s~^ for the 1.173 MeV 
line integrated over the central radian of Galactic longitude. Subsequent analysis of this 
result suggested an upper limit to the mass of 1.7 Mq of '^°Fe in the interstellar medium 
today, which is close to what stellar-chemical evolution calculations yield (Timmes et al. 
1995b). 



4.1 «°Co in SN 1987A 



In old supernova remnants the decay of '^°Co occurs at the same rate as the '^°Fe 
decay that feeds it while in young remnants, like SN 1987A, '^°Co needs to be considered 
independently. The 7-ray luminosity from the beta decay of ^°Co is 

Leo = 6.7 X 10'^ [1 - exp (-KeoMto/t f) + 0.42] X 

, , , / M(«°Co) \ 
exp(-t/r,o) ergs- , (4) 

where (po = 7.0 X lO'* g cm~^ is the column depth at some fiducial time to = 10'^ s, and 
Keo = 0.04 cm^ g~^ is an estimate of the effective opacity (e.g., Woosley, Pinto & Hartmann 
1989). 

The 7-ray deposition efficiency (first term in equation 4) is a key uncertainty. Both 
the column depth and effective opacity are uncertain due to mixing and clumping. A 
more rigorous light curve calculation would distribute the '^°Co, perform Monte Carlo 
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calculations and average over scattering angles (Woosley, Pinto & Hartmann 1989, Woosley 
& Eastman 1994). In the simple approach used here, we absorb all these uncertainties into 
the value for Hqq and examine the sensitivity of the main results given below by varying 
the effective opacity by a factor of two (in addition to the uncertainty associated with the 
yield of «°Co itself). 

Figs. 3a and 3b show the effects expected on the bolometric light curve of SN 1987A 
for an ejected mass of MC^^Co) = 2.0 X 10~® Mq. The figures suggest that '^°Co could have 
been contributing appreciably to the light curve at 1500 days and might even contribute 
to the light curve at the 10% level at 3500 days, but the level of the effect is less than 
the substantial error bars on the current luminosity measurements. There may also be 
significant modifications to the late time light curve due to the freezing out of atomic 
processes and the failure of the steady state assumption (Clayton et al. 1992; Fransson 
& Kozma 1993). The effect of varying the canonical value of the effective opacity (dotted 
curve), /€44=/€6o=0.40 cm^ g~\ by factors of two (dashed curves) is shown in the upper 
panel of Figure 4. The changes are larger for '^°Co than '^'^Ti since '^°Co deposits a larger 
energy per decay than ■*'*Ti. 

For a distance of 50 kpc to SN 1987A and a half-life of 5.27 yr, an ejected mass of 
MC^^Co) = 2 X 10~® would produce a 7-ray line flux of ~ 3 - 5 X 10"'^ photons cm~^ s~^ 
at 1500 days. This line flux might be detectable by future gamma-ray astronomy missions, 
such as INTEGRAL, though the short half-life makes a detection of SN 1987A unlikely. 

5. Summary 

Recent grids of Type II and Type lb models predict '^'^Ti yields of between 10~® to 
lO"'* Mq, with typical values ~ 3 X 10"® Mq for the Type II models and ~ 6 X 10"® Mq 
for the Type lb models. Producing more than lO"'* Mq of '^'^Ti is difficult but making 
much less is easy since '^'^Ti is synthesized in a region very close to the mass-cut. Yields 
of '^°Co from core collapse supernovae vary between 10"'^ and lO"'* Mq, with typical values 
~ 3 X 10"® Mq for the Type II models and ~ 2 X 10"® Mq for the Type lb models. 
Detailed chemical evolution models that use only Type II and lb supernovae and Chan- 
drasekhar mass models for Type la underproduce the solar abundance of ■*'*Ca by a factor 
of 3. The production site of this isotope may be the explosion of sub-Chandrasekhar mass 
white dwarfs which have accreted a helium shell. There is the possibility that a signiflcant 
fraction, if not most, of the steady state diffuse component of the 511 keV emission is from 
the annihilation of positrons due to the decay of ■*'*Ti. 

The line flux derived from COMPTEL and OSSE observations of Cas A give an upper 
bound for the mass of '^'^Ti ejected which is consistent with the core collapse models. The 
abundance of '^'^Ti suggested by the COMPTEL observations raises the possibility that Cas 
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A may have been a brighter supernova than previously thought. If Cas A was a Type lb 
supernovae which experienced significant mass loss that condensed into dust grains prior 
to the explosion, the visual extinction could have been very large. Our estimates suggest 
that an extinction of 10 magnitudes is necessary to understand the historical limits on 
brightness. The abundances of '^'^Ti and '^°Co predicted by Type II supernova models 
suggest the mass of '^'^Ti ejected by SN 1987A may be detectable by the next generation of 
7-ray telescopes, and that '^°Co can contribute to the late time light curves of SN 1987A 
and other Type II supernovae at the 10% level. 
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TABLE 1 

Massive Star ■*'*Ti Yields 



Mass 


±-JJ\. ±±Cai± LjW J. lOLiWll 


1 ' 1 Lj CVJ. 


11 


6.2E-5 


6.2E-5 


12 


8.3E-5 


8.3E-5 


13 


7.1E-5 


6.3E-5 


15 


9.1E-5 


5.7E-5 


18 


8.3E-5 


8.5E-6 


19 


1.2E-4 


1.9E-5 


20 


1.2E-4 


1.4E-5 


22 


1.5E-4 


6.1E-5 


25 


1.6E-4 


3.0E-5 


30 


1.7E-4 


1.4E-4 


35 


1.9E-4 


1.8E-9 


40 


2.0E-4 


1.5E-10 



All entries in solar masses for solar metallicity stars. 
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Figures and Captions 
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1. — Yields of ■*'*Ti and '^°Co from the solar metallicity (a) Type II supernovae models 
of Woosley & Weaver (1995) and (b) Type lb supernovae models of Woosley, Langer 
& Weaver (1995). The large non-monotonic variation in the '^'^Ti yields for the Type II 
supernovae shows the sensitivity to "fall back" for this species made in the innermost 
shell to be ejected. Synthesis of '^°Co is less sensitive to this uncertainty. 



24 



-8 

2 3 



4 5 6 

Interior mass 



7 8 



■4 I I I I I I I I I I I I I I I I I I I I I I I I I I I r. I I I I I I 




Interior mass 

Fig. 2. — Abundance profiles of '^°Co and '^°Fe vs interior mass (Mq) for a solar metallicity, 
25 Mq (a) presupernova star and (b) exploded star at 2.5 X 10'* s. 
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3. — Bolometric light curve of SN 1987A from (a) 500 to 1500 days and (b) 500 to 
3500 days. The observed bolometric luminosity is shown as the filled circles (Suntzeff 
et al. 1992). Contributions from the decay of ®'^'®^Co, ^^Na, ■*'*Ti, and '^°Co are shown 
by the labeled dashed lines while the solid line is the total calculated 7-ray luminosity. 
Ejected masses of M{^^Ti) = 1.0 X 10"^ Mq and M(«°Co) = 2.0 x 10"^ Mq, along 
with effective opacities of = 0.04 cm^ g ^ were used. The relatively long 

half-life and local deposition of the '^''Sc to '^''Ca positron kinetic energy, '^'^Ti tends to 
dominate the bolometric light curve after about 1500 days. However, '^°Co might have 
been contributing appreciably to the light curve at 1500 days and may contribute to 
the bolometric light curve at about the 10% level at 3500 days. 
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4. — Bolometric light curve of SN 1987a from 500 to 3500 days and the effect of varying 
the canonical value of the effective opacity (dotted curve), /€44=/€6o=0.40 cm^ g~\ by 
factors of two (dashed curves). The changes are largest for '^°Co since it deposits a 
larger energy per decay than ■*'*Ti. 
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